MIT-CTP #3347 



Reversible simulation of bipartite product Hamiltonians 



o 
o 

(N 



> 

&\ 
O 

m 
o 
m 
o 

S3 
+-> • 

=3 

CT 1 - 



13 



Andrew M. Childs , X 'E Debbie W. Leung, 2 -Q and Guifre Vidal 2 -© 

1 Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
e Institute for Quantum Information, California Institute of Technology, Pasadena, CA 91125, USA 

(14 March 2003) 

Consider two quantum systems A and B interacting according to a product Hamiltonian H = 
Ha®Hb- We show that any two such Hamiltonians can be used to simulate each other reversibly (i.e., 
without efficiency losses) with the help of local unitary operations and local ancillas. Accordingly, all 
non-local features of a product Hamiltonian — including the rate at which it can be used to produce 
entanglement, transmit classical or quantum information, or simulate other Hamiltonians — depend 
only upon a single parameter. We identify this parameter and use it to obtain an explicit expression 
for the entanglement capacity of all product Hamiltonians. Finally, we show how the notion of 
simulation leads to a natural formulation of measures of the strength of a nonlocal Hamiltonian. 



I. INTRODUCTION 

Suppose two quantum systems A and B are coupled by 
some nontrivial interaction Hamiltonian H ^ Ha ®Ib + 
Ia^Hb- Such a Hamiltonian can be used for a variety of 
information processing tasks, such as transmitting clas- 
sical or quantum information, creating entanglement, or 
simulating other nonlocal evolutions. One of the goals of 
quantum information theory is to quantify the capacity 
of an interaction to perform such information processing 
tasks. 

The various nonlocal properties of Hamiltonians can be 
analyzed in different ways, but they are typically studied 
under a common framework of perfect local control. To 
focus on the uniquely nonlocal features of Hamiltonians, 
all local abilities are regarded as free resources. This 
includes the use of local ancillary degrees of freedom and 
the ability to perform fast local operations to modify the 
evolution. 

One of the tasks mentioned above consists of using H 
and local operations to simulate another Hamiltonian H' , 
a control technique that allows one to modify a naturally 
available interaction into a more convenient one. Any 
bipartite interaction Hamiltonian can simulate any other 
at some nonzero rate 

HUSH However, we would ul- 
timately like to know the most efficient way to use H to 
simulate H' , as well as the optimal rate 1h'\H at which 
the simulation can be accomplished. A method for op- 
timal simulation of two-qubit Hamiltonians is given in 
, and the optimal rate in this case can be expressed in 
terms of a majorization condition However, to the 
best of our knowledge, no optimal simulation rates have 
been reported beyond the two-qubit case. 

Understanding Hamiltonian simulation also provides 
insight into capacities for other information processing 
tasks. Let Ch denote the capacity of Hamiltonian H 
to accomplish one of these tasks, again assuming perfect 
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local control. If Hamiltonian H can be used to simulate 
H' at a rate 1h'\h, then 



Ch > 1h'\h Ch 1 



(1) 



since one could first use H to simulate H' and then use 
H' to accomplish the task. Equation Q is a lower bound 
on the capacity of H, or equivalently, an upper bound on 
the capacity of H' . Of course, such bounds need not 
be tight. For example, the majorization condition for 
optimal simulation of two-qubit Hamiltonians only pro- 
vides a partial order on these Hamiltonians, and thus 
the resulting bounds on capacities — for example, on the 
entanglement capacity [USB — are not always tight. 

The fact that any nonlocal Hamiltonian can simulate 
any other at some nonzero rate means that all interac- 
tions are qualitatively equivalent. A much stronger, quan- 
titative notion of equivalence between interactions comes 
from the possibility of performing a reversible simulation. 
We say that H and W can simulate each other reversibly 
if we can use H to simulate H', and then use H' to sim- 
ulate H back, with no overall loss in efficiency. In terms 
of simulation rates, reversible simulation amounts to the 
condition 

1h\h>1h>\h = 1. (2) 

Notice that if two Hamiltonians H and H' can simulate 
each other reversibly, then their capacities are related by 



Ch — 1h'\h Ch> 



(3) 



as can be seen by applying (JIJ in both directions. Fur- 
thermore, if every pair of Hamiltonians in some given set 
can simulate each other reversibly, then simulation pro- 
vides a total order on the set. Thus the nonlocal proper- 
ties of the entire set can be studied by focusing on only 
one Hamiltonian in the set. 

In this paper we consider the set of bipartite Hamil- 
tonians that can be written as a tensor product of the 
form 



H = H A ®H £ 



(4) 



where Ha acts on system A and Hb acts on system B. 
We shall call such a Hamiltonian a product Hamiltonian 



2 



for short. An example of a product Hamiltonian in a 
two-qubit system is the Ising interaction Hi s i n „, 



1 

-1 



(5) 



Our main result is an explicit protocol for the reversible 
simulation of any product Hamiltonian by another. It 
follows that the nonlocal properties of a product Hamil- 
tonian H depend entirely on a single parameter. We 
denote this parameter by K®(H), and choose it to be 
the rate 7ij Isin \h a t which H can simulate the Ising in- 
teraction. We find that 



(0) 



where (As) denotes the difference between the 

largest and the smallest eigenvalues of Ha (Hb)- The 
optimal simulation rate between any two product Hamil- 
tonians H and H' can be written in terms of Ka as 



1h'\h 



(7) 



so that any capacity Ch known for just one product 
Hamiltonian can be easily computed for any other prod- 
uct Hamiltonian using @ and Q. In particular, we use 
previous results for the Ising interaction to obtain a 
simple expression for the entanglement capacity of any 
product Hamiltonian. 

In addition to quantifying the ability of an interac- 
tion to perform particular tasks, one can imagine defin- 
ing abstract measures of the nonlocality of an interac- 
tion. Reference introduced the notion of a strength 
measure K(U) of a nonlocal quantum operation. Three 
axioms that every strength measure should satisfy were 
proposed, along with several additional desirable proper- 
ties. Here we consider strength measures K(H) for non- 
local Hamiltonians. We will formulate a single axiom in 
terms of Hamiltonian simulation that implies many de- 
sirable properties. For the special case of product Hamil- 
tonians, our results imply an essentially unique measure, 
K®(H) in JBJ, with several additional properties. 

The structure of the paper is as follows. In Section ITU 
we describe the problem of bipartite Hamiltonian simu- 
lation in more detail and discuss a set of basic simulation 
rules that can be used to build up all possible simula- 
tions by composition. In Section HTT1 we derive our main 
result, namely that all tensor product Hamiltonians can 
reversibly simulate each other. In Section llVI we apply 
this result to the calculation of entanglement capacities, 
and in Section [V] we relate Hamiltonian simulation to 
the strength measure formalism. Finally, in Section IVII 
we conclude with a discussion of our findings and of some 
open problems. 



II. SIMULATING BIPARTITE HAMILTONIANS 

The problem of bipartite Hamiltonian simulation can 
be posed as follows. We consider two quantum systems 



A and B held by Alice and Bob, respectively. The sys- 
tems interact according to some nonlocal Hamiltonian 
H 7^ Ha ®Ib+Ia®Hb- Alice and Bob want to use H to 
produce an evolution according to some other bipartite 
Hamiltonian H' . In order to do so, they are given the 
additional resources of synchronized clocks and perfect 
local control. They may attach or discard local ancillary 
systems and they may apply (arbitrarily fast) local opera- 
tions, which can be assumed to be unitary without loss of 
generality 0- As a side remark, classical communication 
between Alice and Bob would not increase the optimal 
simulation rate . Note that the goal of the simulation 
is not to produce the Hamiltonian evolution e~ lH * for 
a particular time t, but rather to stroboscopically track 
the evolution e~ lH * for arbitrarily closely spaced values 
of time. A detailed formulation of the problem can be 
found in 0. 

We next present a list of rules for nonlocal Hamilto- 
nian simulation. By composition, these rules give rise to 
all possible simulations achievable with local operations 
and ancillary systems. We present five basic rules, as 
well as three additional rules that can be obtained by 
combining the basic ones. We use the shorthand nota- 
tion H — > H' to represent the possibility of simulating 
H' by H at the rate 7h'\h, and the notation H < — > H' 
to indicate that, in addition, the simulation can be re- 
versed without overall efficiency losses, as in (|2J). We say 
that two Hamiltonians are locally equivalent if they can 
simulate each other reversibly at unit rate. 

The first two basic rules merely make precise the notion 
of Hamiltonian evolution. They do not involve any oper- 
ational procedure, nor assume any ability to control the 
system. The first rule makes precise the notion of rescal- 
ing the evolution time: a Hamiltonian H can reversibly 
simulate another Hamiltonian H' = cH that only differs 
by a positive multiplicative constant c. 

Rule 1 (Rescaling). For any c > 0, 



H 



cH . 



7 



c.H\H 



(8) 



Note that it is important that c > 0. In general, Hamil- 
tonians H and —H cannot simulate each other reversibly 
(see 00 for examples). 

The second rule makes precise what it means for a 
Hamiltonian to act on a subsystem. In the bipartite 
setting, the complete system can be described by sub- 
systems A, B on which H acts and ancillary subsystems 
A', B' on which it acts trivially. 

Rule 2 (Ancillas). For any dimension of the ancillary 
Hilbert space H.a> <£> T~Cb' , 



H 



H <g> I A 'B' 



,\H 



1. 



(9) 



The next two basic rules arise from the possibility of 
switching on local Hamiltonians. 

Rule 3 (Local Hamiltonians). Any local Hamiltonian 
of the form Hq = Ha (^Ib+Ia® Hb can be produced at 
no cost. 
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Recall that for any unitary transformation U, we have 
U e -iHtjj\ _ e -iUHU t_ Therefore, by means of local 
unitaries, any Hamiltonian H is locally equivalent to any 
other that is obtained from it by local unitary conjuga- 
tion. 

Rule 4 (Local unitaries). For any local unitary oper- 
ation U = Ua®Ub, 



H < — ► UHU^ . 



1uhW\h 



1 . 



(10) 



Rules allow Alice and Bob to produce Hamilto- 
nians that differ from the original interaction H. The 
Lie- Trotter product formula, 



tell us how two of these Hamiltonians can be combined 
into a new one by alternately simulating each of them 
individually. With the help of Rule ^ we obtain the last 
basic rule. 

Rule 5 (Convex combination). For any H\,H 2 and 
< p < 1 } the simulation 



(l-p)H 2 



H' = pH t + (1 - p)H 2 



(12) 



is possible with rate 1h'\pH 1 ;(i-p)h 2 ^ 1- 



Here we have considered the use of Hamiltonian H\ for 
a fraction of time p and Hamiltonian H 2 for a fraction 
of time 1 — p, and the rate of simulating H' is computed 
by adding these two times together. Let us stress that 
(|12f> assumes only the local ability to switch on and off 
the constituent Hamiltonians, and that only one Hamil- 
tonian is acting at a time. Notice that this is the only 
basic simulation rule where irreversibility may occur. Al- 
though Alice and Bob can use H' to simulate back Hi 
and H 2 , in general they will incur an overall loss in effi- 
ciency by doing so. 

These basic rules can be combined in various ways. We 
state three particularly useful combinations as additional 
rules. First, from Rules [3] and [5] a local part Hq can be 
added to the given nonlocal Hamiltonian reversibly. 



Rule 6 (Adding a local Hamiltonian). 



H 



H + Hn 
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H+H \H 



= 1 . 



(13) 



Second, local unitary conjugation and convex combi- 
nation can be composed into what we shall call a local 
unitary mixing of H. 

Rule 7 (Local unitary mixing). For any set of lo- 
cal unitary transformations Ui = Uas ® Ubs and any 
probability distribution pi (pi > and J2iPi = -U.- 



H 



1h>\h>1- (14) 



Note that Rules [3][7| are stated without assuming local 
control over the ancillas. In the two-qubit case, Rules 
n EH and [7| describe all relevant simulations because lo- 
cal control over ancillas is known to be unnecessary for 
optimal simulations |2|. 

More generally, we allow local control over ancillas in 
our simulation model. By Rule [3 Rules |3J{7| can be ex- 
tended to include ancillas as well. Control over ancillas 
gives extra freedom in the simulation, and is known to 
improve the achievable simulation rates in some cases p| . 

Our last rule is concerned with any simulation in which 
the original Hamiltonian H and the simulated Hamilto- 
nian H' act on systems with different dimensions. Let 
Ti-A <8> Hb and H.a> <8> Ti-B' denote the Hilbert spaces 
on which H and H' act, with dimensions dA-,dB and 
dA'yds 1 , where dA > and ds > ds>- For simplicity, 
we assume that H = Ha <8> Hb is a product Hamiltonian. 
If it were not, then we could expand H as a linear com- 
bination of product Hamiltonians, H = J2i Ha,% <S> Hb,i, 
and the following would hold for each of the terms in 
the expansion. Let vectors \j)a (1 < j < c?a) denote an 
orthonormal basis in Ha- We can express Ha as 



H A = 



J,, ct 
C Jj 



(15) 



where J\\ is the restriction of Ha onto the subspace 
— s P anne d by vectors \j)a (1 < j < dA'), and 
J_l the projection onto its orthogonal complement. Con- 
sider also an analogous decomposition for Hb- Then we 
have the following: 

Rule 8 (Reduction to a local subspace). The sim- 
ulation 



H = 



' J \\ 


0' 






c 






D K x 



H' = J\\®K\\ (16) 



is possible with rate 1h'\h ^ 1- 

The last rule can be obtained by using the following 
lemma twice. The lemma shows how to simulate J|| ®Hb 
using Ha ® H B ■ 



Lemma 1 

H = 



The simulation 
J11 C f 



H 



B 



H' = J,, $ H E 



(17) 



C Ji 

is possible with rate 1h'\H ^ !• 

Proof. We divide the simulation into two steps, (i) First, 
by unitary mixing (Rule[7|) with 





1 








Pi = 


2' 




= Ia®i b ; 




1 


u 2 




"In 


P2 = 


2' 




-i± 



If 



(18) 



where /|| and I± denote restrictions of the identity oper- 
ator, we achieve the simulation 



H = 



J|| ct 
C Ji 



H f 



H" = 



J,, 

Ji 



H b (19) 
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with unit rate, so 1h"\h > 1- (H) Second, we use H" to 
simulate H' = Jii i?s as follows. Suppose the goal is 
to evolve \ip}A'B according to e~ lH * . We assume system 
A is in state by local control. Therefore the joint 
state of systems AA'B is initially \1) A \ip) A , B . Let Vaa' 
denote a unitary transformation such that 

Vaa>\1)a\3)a> = \3)a\ 1 )a" l<3<d A >. (20) 

Then the following three steps can be used to complete 
the desired simulation: 

1. We apply unitary Vaa'i placing \ip) A , B in the sub- 
space Ha\\ <8> Kb C Ha ® Hb- 

2. We make AB evolve according to H". Notice that 
at all times e~ %H t \ip) AB is supported in Ha\\ ®Wb, 
and that H" acts on this subspace as J\\ <E> Hb- 

3. We reverse Vaa 1 , so that the net evolution on A'B 
has been e~ lH t \ip) A , B - 

This completes the proof. □ 

III. REVERSIBLE SIMULATION 

In this section, we present the main result of the paper, 
the reversible simulation of tensor product Hamiltonians 
H = Ha <8> Hb- We will consider product Hamiltoni- 
ans in a certain standard form. Using Rule we may 
diagonalize Ha and Hb, so we need only consider their 
eigenvalues. It will also be convenient to modify Ha so 
that the largest and smallest eigenvalues, A™ ax and A™ 1 ", 
are equal in magnitude, and similarly for Hb- This can 
be done by adding a term proportional to the identity to 
each of Ha and Hb, i.e., 

(H A + cl) <gs (H B + dl) = H A ® H B + cl <g> H B 

+ dH A ® I + cdl (g) I . (21) 

The resulting Hamiltonian is locally equivalent to H since 
they differ only by local terms (Rule[SJ). Furthermore, 
since 

(cH A ) ® (H B /c) = H a ®H b , (22) 

we may assume A^ ax - A^ in = A"g ax - Ag in = A without 
loss of generality. 

Having put all product Hamiltonians into a standard 
form, we are ready to show that they can reversibly sim- 
ulate each other. By the transitivity of reversible simula- 
tion, it suffices to show that all product Hamiltonians can 
reversibly simulate the Ising interaction -ffising = o~ z ®<j z . 

Lemma 2. Any tensor product Hamiltonian H = Ha <S> 
Hb can reversibly simulate the Ising interaction 

H^H lsine , l Hlsins]H = Ja 2 . (23) 



Proof. For any nonlocal H\ and Hi, we have 

'yH 1 \H^H 2 \H 1 < 1. (24) 

Otherwise we could use Hi to simulate itself with sim- 
ulation rate greater than 1, which is a contradiction. 1 
It thus suffices to show that 7ff Ising |_ff > A 2 /4 and 

l H \ HlBins > 4/A 2 . 

Let H act on a Hilbert space TLa®Hb with dimensions 
dA and ds- Since H = Ha ® Hb is in the standard form, 
we may write 

H = |A 2 diag(l, a 2 , a 3 , . . . , a dA -i, -1) 

®diag(l,& 2l &3,... ) b<fe-i,-l) (25) 

where 

1 = ai > <i2 > ... > aa A = — 1 (26) 
1 = h > b 2 > ... >b dB = -l, (27) 

and the corresponding eigenvectors are \ j) A (1 < j ' < c?a) 
and \j) B (1 < j < d B ). 

We can simulate the Ising interaction using H by re- 
stricting to the subspace spanned by the extremal eigen- 
vectors of Ha and H B (\1)a\1)b, \d A )A\l)B, 
\dA)A\dB)B) according to Rule |H| In this subspace, H 
acts as (A 2 /4)ifisi ng . Therefore we have 

1h 1s ^\h> A 2 /4. (28) 

In order to show how to use the Ising interaction -ffising 
to simulate H, we consider a concatenation of two simu- 
lations, 

-falsing — ► H" — ► H . (29) 

Here H" = jA 2 H' AA ,®)H BB , acts on local Hilbert spaces 
of dimensions 2d a and 2ds , and reads 

H" = |A 2 diag(l, -1, as, -a 2 , •••,-1,1) 

®diag(l, -1,62,-62,.-., -1,1). (30) 

Clearly, we can use Rule to simulate H by H" with 
unit simulation rate. Therefore we need only focus on 
the simulation of H" by -fusing- In turn, this can be 
decomposed into two similar simulations, 

-falsing — ► H AA , ® a z — > H' AA , ® H BB , , (31) 

each one with unit rate. In order to simulate H' AA , ® o~ z 
using a z ® a z , we append a d^-dimensional ancilla A 



If 1 h^h-z'Y H 2 \H! > 1, then we could concatenate several simula- 
tions Hi > H2 > Hi > • ■ ■ > H'i > Hi to obtain that 

the optimal simulation rate Ih^h-l ls infinite. Recalling that 
any bipartite nonlocal Hamiltonian can simulate any other one 
at finite rate, we would conclude that 1 h\h i s a l so infinite for any 
bipartite Hamiltonian. This would contradict, for instance, the 
results of Q showing that 1h\h = 1 f° r an nonlocal two-qubit 
Hamiltonians. 
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to qubit A! (with H\ s 
Hamiltonian 



acting on A'B') to obtain the 



We define 



so that 



= (I ® <Tz,A>) ® <Tz (32) 

= diag(l,-l,l,-l,...,l,-l)®o*. (33) 

Pj = (aj + l)/2 (34) 

1 = Pi > P2 > ■ ■ ■ > Pd A = . (35) 



Furthermore, we define d,A local unitary operations 
where U 3 exchanges the (2j — l)-th and (2j)- 
th basis vectors of AA! . To evolve under H' AA , ® cr z for 
a small time S, we apply each Uj at time t = pj5 and [/"• 
at time t = 5. Equivalcntly, we can use Rule □ with an 
appropriate probability distribution and set of unitaries. 
Thus we can use ffising to simulate H AA , (g> a z with unit 
efficiency. The second simulation in <|31[) is achieved sim- 
ilarly. The overall rate for -ffising to simulate ff" or ff is 
thus 4/ A 2 by Ruled □ 

We have shown that any product Hamiltonian H can 
reversibly simulate the Ising interaction ffising with rate 
7ff Islng |ff = K®{H), where 



KafH) 



(36) 



Our main result then follows: 



Theorem 3. Any product Hamiltonian H can reversibly 
simulate any other product Hamiltonian H' , with simu- 
lation rate given by 



~<H'\H 



(37) 



As discussed previously, in general a bipartite Hamilto- 
nian H cannot reversibly simulate — H . Similarly in gen- 
eral H cannot reversibly simulate its complex conjugate 
H*, nor the Hamiltonian ff° resulting from swapping 
systems A and B. However, for product Hamiltonians, 
all these Hamiltonians are locally equivalent. 

Corollary 4. For any product Hamiltonian H , 

K®(H) = K 9 (-H) = K®(H*) = K®{H«) . (38) 

Theorem |31 can be extended to the case of a sum of bi- 
partite product Hamiltonians acting on separate systems. 
If Hi and H2 are two Hamiltonians acting, respectively, 
on bipartite systems A\B\ and A2B2, we let Hi EE H2 
denote their sum. 2 In fact, Hi EH H2 can reversibly simu- 
late a product Hamiltonian H acting on a single bipartite 
system AB. 



2 We use the symbol EB rather than + to emphasize that the Hamil- 
tonians being summed act on different pairs of systems. In other 
words, H AB E H A , B , = H AB ® I A , B , + I AB ® H A , B , . 



Lemma 5. If Hi, H2, and H' are product Hamiltonians, 
the simulation 

Hi EH H 2 < — ► H' (39) 
can be achieved reversibly, with simulation rate 

^H'\ Hl \SH 2 = lw\ Hl + ~Ih>\h 2 • (40) 

Proof. Because of Theorem |3 we only need to show that 
the Hamiltonian H' = cHi s - mg EH d -ffising, c,d € M, can 
reversibly simulate H = (|c| + |d|) ffising at unit rate. In 
addition, Corollary 0] implies that we need only consider 
the case c, d > 0. 

By Rule |5J ff = (c + rf)ffising is locally equivalent to 



Ji = (c + d) a zA (g) a zB ® I A'B' 



(41) 



In turn, using local unitaries to swap A with A! and B 
with B' (Rule^J), Ji is locally equivalent to 



■h 



d) I A! 



a zA ' ® o- zB > 



(42) 



Then we can simulate ff' = [c/(c+d)] Ji + [d/(c+d)) J 2 by 
convex combination (Rule|SJ) of Ji and J2, which shows 
that 1 H <\H > !• 

For the reverse simulation, note that ff' is locally 
equivalent to each of the following four Hamiltonians: 

c(a zA <£> I A' ®o- zB ® Ib> ) + d(I A ®o- zA '® Ib ®e zB ') , 
c( I a ®o- zA '®cr zB ® I B , ) + d(a zA ® I A' ® Ib ®ct z b>) , 
c(a zA ® Ia> ® Ib ®ct zB ') + d( I A ®(T Z A'®cr zB ® I B > ) , 
c( I a ®<J Z A'® Ib ®VzB>) + d(a z A® I A' ®u zB ® I B ' ) ■ 

(43) 

Each of these Hamiltonians can be obtained from ff ' ac- 
cording to Rule 0] by swapping A with A' and B with B' 
as necessary. An equally weighted convex combination of 
these four Hamiltonians gives, after rearranging terms, 

c+d 

— — {<j zA ®Ia' + Ia®cta> ) ® (o- z b®Ib> + Ib®o-b< ) , (44) 



a tensor product Hamiltonian with A 2 = 4 (c+d). There- 
fore 1h\h' > 1, which completes the proof. □ 

It follows from Lemma[5]that is additive under the 
sum of product Hamiltonians acting on different pairs of 
systems, 



K®(Hi EH ff 2 ) = K 9 (Hx) + K®{H 2 ) . 



(45) 



More generally, for ff = EBj ff; and ff' = EBjff 4 ', where 
all Hi and ff • are bipartite product Hamiltonians, we can 
perform the simulation 



ff = fflff. 



ff ' = ffl Hi 



reversibly, with simulation rate given by 



1h>\h 



(46) 



(47) 
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Finally, in the Appendix, we present a case of reversible 
Hamiltonian simulation that is possible when in addition 
to local operations and ancillas, catalytic pre-shared en- 
tanglement is available to Alice and Bob. The simulation 
can be made reversible only in the presence of entangle- 
ment, but the entanglement is not used up during the 
simulation 



IV. ENTANGLEMENT CAPACITY 

In this section we use the results on reversible simula- 
tion from Sec. lIIII to determine the optimal way a product 
Hamiltonian can be used to produce entanglement. 

The problem of optimal entanglement generation by an 
interaction Hamiltonian has been approached in different 
ways. In |(| the single-shot scenario was considered. This 
corresponds to the setting in which an interaction H is 
used, with the help of fast local unitary transformations 
(but without control over local ancillas), to maximize the 
rate of increase of entanglement between a single copy 
of the two interacting systems. This situation is of in- 
terest for many present day experiments aiming to pro- 
duce entangled states in quantum optics, nuclear mag- 
netic resonance, or condensed matter systems ^l|- The 
optimal rate at which entanglement can be generated in 
the single-shot scenario is known as the entanglement ca- 
pability of H, Th, and its value has been determined for 
any two-qubit Hamiltonian 6] . For the Ising interaction 
it reads 



Hi. 



= a w 1.9123. 



where 



2 max \/x(l 

xe[o,i] 



a;) log 



1-x 



(48) 



(49) 



with the maximum obtained at xq « 0.9168. 

In contrast, Q considers entanglement generation in 
the asymptotic scenario, where many copies of the in- 
teracting systems (and local ancillas) can be used col- 
lectively to produce entanglement, possibly at a higher 
rate than in the single-shot case. The asymptotic entan- 
glement capacity (or entanglement capacity for short) of 
interaction H , denoted Eh, is of interest in the context of 
understanding the ultimate limitations of quantum me- 
chanical systems to process quantum information. Ref- 
erence showed that entanglement capacities equal en- 
tanglement capabilities when ancillas are allowed in the 
single-shot setting. 

Finally, it was shown in Q that for a number of two- 
qubit Hamiltonians, all single-shot and asymptotic sce- 
narios lead to the same optimal entanglement generation 
rates. In particular, 



E 



Hi, 



Hz, 



a . 



(50) 



Combining Q and Theorem |3 we obtain an expres- 
sion for the entanglement capacity of any product Hamil- 
tonian H, since we have Eh = K^(H) Ejj lsin ■ 



Theorem 6. For any product Hamiltonian H , 



E H = -A A A B . 



(51) 



Similarly, we can compute E Hl mH 2 f° r the sum of two 
product Hamiltonians H\ and H^ acting on different 
pairs of systems: 



E 



H t mH 2 — + Eh. 



(52) 



In fact, H51fl also corresponds to the single-shot capa- 
bility Th, since it can be obtained without using ancillas. 



Corollary 7. For any product Hamiltonian H , 

T H = ^A A A B . 
Proof. The explicit optimal input state is 



(53) 



Vxo\+)a <8> \+) B + Wl - x \-)a ® \-)b (54) 



where 



\±) A = -^(\l) A ±\d A ) A ) 



(55) 



and similarly for system B. Here |1)^ and \d A ) A repre- 
sent the eigenstates of H A corresponding to the largest 
and smallest eigenvalues, respectively. That this state 
achieves Tjj can be seen by substitution into Eqs. (17) 
and (18) of 0. □ 

Likewise, Q52JI can also be achieved without ancillas, be- 
cause the protocol used in Lemma |3] does not involve 
ancillas: 



FhiBH 2 = Tffi + Th 2 



(56) 



As a side remark, [l2j has reported the restricted case 
of 1511) in which H A and Hb have eigenvalues ±1, but 
this property has no special significance; only the tensor 
product structure is important. 

Equation (|51|l could also be proved directly using es- 
sentially the same arguments that appear in Using 
the observation that the matrix with elements 



A\(H A ) jk (H B ) jk \/A 2 



(57) 



(where the matrix elements of H A and Hb are taken 
in any orthonormal basis) is doubly substochastic, an 
upper bound on E H can be obtained just as in Eq. (26) of 
[8J. However, the approach using reversible simulation is 
more powerful, since it applies to any capacity, even those 
for which the capacity of Hi s - mg is yet to be found, such 
as the capacity for communicating classical information. 

Finally, we note that Theorem [5] can be extended to 
Hamiltonians that can be reversibly simulated using cat- 
alytic entanglement, such as those mentioned in Sec. IIIII 
and further described in the Appendix. This class of 
Hamiltonians includes, as a special case, the full set of 
two-qubit Hamiltonians of the form [i x a x ®a x -\-iJLy o~ v ®o~ v 
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considered in |8( ■ In the context of asymptotic entangle- 
ment capacity, catalytic resources need not be considered 
as additional requirements since the cost of first obtaining 
any catalytic resource can be made negligible How- 
ever, it turns out that for the Hamiltonians discussed 
in the Appendix, catalytic entanglement is actually not 
necessary to achieve the entanglement capacity. 

V. STRENGTH MEASURES FOR BIPARTITE 
HAMILTONIANS 

In this section, we show that the notion of nonlocal 
Hamiltonian simulation can be used to define measures 
of the strength of a nonlocal Hamiltonian. We proceed 
along the lines of 0, which introduces a formalism of 
strength measures for nonlocal quantum operations, with 
the goal of quantifying their nonlocality. Three neces- 
sary axioms and a number of other desirable properties 
for such measures were proposed. Here we show that in 
the case of nonlocal Hamiltonians, a single axiom implies 
many of these properties. 

We denote a strength measure for Hamiltonian H as 
K{H). The only requirement we impose on K(H) is 
that it does not increase under Hamiltonian simulation 
by local manipulations. That is, if H can be used to 
simulate H' at a rate 1h'\Hi then the strength measure 
for 1 h'\h H' should be no greater than that for H . This 
is motivated by the idea that any measure of the non- 
locality of an interaction should not increase under local 
manipulations. 

Axiom (Monotonicity). Any strength measure K(H) 
must satisfy 

K(H) > 1 H , W K(H') (58) 
for any two Hamiltonians H,H'. 3 

In the following, we exclude the trivial strength mea- 
sure K(H) — for all H. Likewise, we exclude any un- 
physical strength measure that is infinite for a bounded 
H. 



3 One might also be interested in defining a strength measure in 
a setting where more than just a single Hamiltonian is avail- 
able. Suppose we can use two switchable Hamiltonians Hi 
and H2 to achieve some task, and that we can only have one 
Hamiltonian switched on at a time. Let p and 1 — p (0 < 
p < 1) characterize the relative frequency with which we use 
each of the two Hamiltonians. Then we can define a func- 
tion K(pHi; (1 — p)H2) as a strength measure. More gener- 
ally, K may have any number of arguments. Suppose we can 
use n Hamiltonians {piHi; P2H2', • • • \PnH n }y where pi (pi > 0, 
EiPi = 1) indicates the relative frequency of Hi, to simu- 
late Hamiltonians {qiH[; qiH'^\ ■ • ■ ; q m H' m }. Let ^ { qi H'}\{ Pi H i } 
denote the rate at which this can be achieved. Then we 
would require that a strength measure K fulfills K({piHi}) > 
{q i H'}\{p i H i }K{{(liH' i }). Several additional properties can be 
derived for K. In particular, Rule l5l implies that K({piHi}) > 



We have already encountered two examples of non- 
trivial functions that do not increase under local manipu- 
lations: (i) any capacity Ch satisfies the Axiom because 
of and (ii) for any fixed target Hamiltonian -fftarget, 
the simulation rate 7m ah satisfies the Axiom as a 

-'-•target 1 11 

function of H . In particular, the function 

K®(H):=1 Huing{H (59) 

used in previous sections is a strength measure. 

A function K(H) satisfying the Axiom automatically 
has a number of properties that we describe below for 
the bipartite case (although much of the discussion can 
be generalized to more than two systems). Several of 
these properties appeared in the original formulation in 
P 

Property 1 (Positivity). 

K{H) > (60) 
with equality if and only if H is local. 
Property 2 (Homogeneity). For any c > 0, 

K{cH) = cK{H) . (61) 

Property 3 (Stability under addition of ancillas). 

For any ancillary system A' B' , 

K(H)=K(H®I A , B ,). (62) 

Property 4 (Local unitary invariance). For any lo- 
cal unitary operation U = Ua®Ur, 

K{H) = K{UHU f ) . (63) 

Property 5 (Invariance under local addition). For 

any local Hamiltonian Hq, 

K(H) = K(H + H Q ) . (64) 

Property 6 (Local unitary mixing). For any set of 

local unitary operators Ui — Uas ® Ub,i and probability 
distribution pi, 

K{H) > K (^T Pl U t HU^ . (65) 

Property 7 (Reduction to a local subspace). Let Jn 

denote the m x m upper left submatrix of J, and Ku the 
m! x m! upper left submatrix of K. Let H — Ji ® Ki 
be a product expansion of a bipartite Hamiltonian H , and 
let H' = Jjii ® Km be its restriction to the upper left 
subspace. Then 

K(H) > K(H') . (66) 

Property 8 (Continuity). For any nonlocal Hamilto- 
nian H and any bounded Hermitian operator J, 

lim K {H + e J) - K(H) = 0. (67) 
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We now derive Properties^S]from the Axiom and our 
simulation rules in Section^ First we prove Property^ 

Proof. Let Hq be any local Hamiltonian and H be any 
bounded nonlocal Hamiltonian. Hq cannot simulate H, 
so ThIHo = 0- The Axiom then implies K{Hq) > 0. 
Through local control, Hq can be simulated without using 
H, so 1h \h > 0- The Axiom then implies K(H) > 
0. However, 1h \h is unbounded. If K(H) is bounded, 
K(H ) must be 0. Finally, if K(H) = 0, the fact that it 
can simulate all other (bounded) Hamiltonians H 1 with 
nonzero rate implies K(H') = 0, which is excluded. Thus 
K{H) > 0. □ 

Properties [2H3 follow from the simulations described, re- 
spectively, in Rules □ El E EHH1 Finally, the proof of 
Property [S] is based on the fact that any two nonlo- 
cal bipartite Hamiltonians can simulate each other with 
nonzero rate: 

Proof. Because nonlocal strength measures are invariant 
under addition of local terms (Property [SJ, we can focus 
on a Hamiltonian H without local terms, that is, tr^ H = 
tiB H = 0, and similarly for J. Consider first the case 
where J is proportional to H. Then Property [21 ensures 
continuity. Suppose now that J is not proportional to H, 
so that H + e J 7^ Ve and is therefore always nonlocal. 
We then consider the following two simulations: 



H 
H 



H 



eJ 



eJ 
• H . 



(68) 
(69) 



A possible procedure for l|68[) is to use H to simulate 
itself at rate 1h\h an d eJ at rate 7 e j|_y, and to add the 
two Hamiltonians together by Rule This gives a lower 
bound on the optimal simulation rate 



^H+c,J\H > 



( 1 



Then, for small e, 



"feJ\H 



= 1 



lj\H 



lj\H 



+ 0{e 2 ). 



(70) 
(71) 



Similarly, for simulation (|69|l we arrive at 



^H\H+tJ > 1 



+ 0(e 2 ). 



-J\H+eJ 



(72) 
) for 



"f-J\H+eJ 

Furthermore, let g := min o<e<eoOjjff > 
any fixed e > 0. Then the Axiom implies 

K{H) > (l-eg )K{H + eJ) + O(e 2 ) (73) 
K{H + eJ) > (l-eg )K{H) + O{e 2 ) (74) 



and 



\K(H) - K(H + eJ)\ <eg K + 0(e 2 ) . 



(75) 



where K := max o<e<e K(H+eJ), which proves the con- 
tinuity oi K. □ 



For bipartite tensor product Hamiltonians, the follow- 
ing additional properties hold. 

Property 9 (Inverse Hamiltonian) . For any product 
Hamiltonian H = Ha ® Hb, 



K{-H) = K(H) 



(76) 



Property 10 (Complex conjugation). For any prod- 
uct Hamiltonian H — Ha ® Hb, 



K(H*) = K(H) . 



(77) 



Property 11 (Exchange). For any product Hamilto- 
nian H = Ha <S> Hb and for H° the same Hamiltonian 
acting on exchanged systems, 



K(H°) = K(H) . 



(78) 



Property 12 (Additivity). For product Hamiltonians 
Hi and H 2 acting on different pairs of systems A\B\ and 
A 2 B 2 , 



K{H X ffl H 2 ) = K [Hi) + K(H 2 ) 



(79) 



Properties l^tllll follow from Corollary 01 Property ITS! 
follows from Lemma 

The above Axiom relates the notion of strength mea- 
sure to that of Hamiltonian simulation. The following 
lemma shows that identifying strength measures is of in- 
terest in order to establish bounds on simulation rates — 
and thereby also bounds on Hamiltonian capacities, as 
discussed in Sec. H] 

Lemma 8. The optimal simulation rate 1h'\h corre- 
sponds to an optimization over all possible strength mea- 
sures K , 



■ K{H) 
7h>\h = nun J^JJjrj 



(80) 



Proof. From the Axiom, 1h'\h < k(h') ^ or an ^ strength 
measure K. In addition, the minimum is achieved by the 
strength measure K(H) = 1h'\h- d 

The many functions that fulfill the Axiom form a con- 
vex cone. In other words, if the functions K\ and K 2 
satisfy the Axiom, then so does c\K\ + c 2 K 2 for any 
ci,c 2 > 0. Of special interest is the subset of extremal 
strength measures, i.e., those that cannot be expressed 
as a positive sum of others, since we can restrict the op- 
timization in H80|) to such functions. Ideally, we would 
like to find a finite subset of extremal strength measures 
that form a complete set, in that optimization over this 
set gives the optimal rate 1h'\h I0r au H> H' . We will de- 
scribe two examples of complete sets of extremal strength 
measures. 

The first case corresponds to the simulation of two- 
qubit Hamiltonians by two-qubit Hamiltonians. Any 
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such Hamiltonian can be written, up to local terms and 
local unitary transformations, as 



H = 



E 

i=x,y,x 



A, 



Az > Aj, > |A 2 



The set of strength measures Ki 



K 3 



A,. 

Arc + Aj, 
Arc + A y + X z , 



A, 



(81) 



(82) 
(83) 
(84) 



is a complete set of extremal strength measures in that 
any optimal rate 1h'\h can be obtained from the opti- 
mization 



1h-\h = 



K^H) K 2 (H) K 3 (H) 



Ki[H'Y K 2 (H'Y K 3 (H>) 



(85) 



The second case corresponds to the simulation of prod- 
uct Hamiltonians by product Hamiltonians. We have 
seen that in this case the simulation can always be made 
reversible. This implies that, up to a multiplicative con- 
stant, there is a unique strength measure for product 
Hamiltonians, 



1 A A 

:= -A A A B 



(86) 



Indeed, because of reversibility, any function K fulfilling 
the Axiom satisfies 



K{H)=1 H , ]H K{H') 
for all product Hamiltonians H,H', so that 
K{H) K®(H) 

— I H'\H — 



(87) 



K{H>) 



VI. DISCUSSION 

We have seen that all tensor product Hamiltonians 
can simulate each other reversibly, so that their nonlo- 
cal properties are characterized entirely by the quantity 
K®(H) given in ©. This is an example of lossless in- 
terconversion of resources, an appealing situation in in- 
formation theory. A related example is the problem of 
communication through a one-way classical channel. By 
Shannon's noisy coding theorem 0] together with the 
reverse Shannon theorem 0], all classical channels can 
simulate each other reversibly (in the presence of free 
shared randomness), and hence they can be character- 
ized entirely in terms of a single quantity, their capac- 
ity. Similarly, in the presence of free shared entangle- 
ment, all one-way quantum channels can simulate each 
other reversibly (at least on certain input ensembles |15|). 
and thus they are characterized entirely in terms of their 
entanglement-assisted capacity for sending classical in- 
formation. 



In Section IIVI we saw how Theorem [3] can be used 
to extend previous results for two-qubit Hamiltonians to 
product Hamiltonians. Another such extension can be 
obtained for the problem of using bipartite Hamiltoni- 
ans to simulate bipartite unitary gates. In the case of 
two-qubit systems, it is known how to optimally pro- 
duce any two-qubit gate using any two-qubit Hamilto- 
nian [To . ITtI Il8| . Since all product Hamiltonians are 
equivalent to some multiple of the Ising interaction, this 
result immediately provides the optimal way to use any 
product Hamiltonian to simulate any two-qubit unitary 
gate, such as the controlled-not gate. 

In view of our results, it will be interesting to improve 
our understanding of the properties of the Ising interac- 
tion. For example, as previously mentioned, a calculation 
of the communication capacity of the Ising interaction 
would provide a formula for the communication capacity 
of all product Hamiltonians. 

It is also interesting to consider Hamiltonian simula- 
tion in the multipartite case. All of the Rules from Sec. ITTI 
have multipartite analogues, and much of the general 
discussion on strength measures for Hamiltonians from 
Sec. [V] can be carried over as well. Furthermore, The- 
orem [3] can be generalized to more than two parties in 
the special case in which the individual tensor factors 
are traceless. However, much remains to be done in the 
general multipartite case. 

Of course, the set of tensor product Hamiltonians is 
clearly a special subset of all bipartite Hamiltonians, and 
thus may not be representative of the general problem 
of bipartite Hamiltonian simulation. For example, we 
have seen that product Hamiltonians admit a total or- 
der, whereas even in the two-qubit case, general Hamil- 
tonians only admit a partial order. Also, note that for 
product Hamiltonians, H and —H are locally equivalent, 
so that in particular, Eh = E^h- However, while this is 
true for all two qubit Hamiltonians, numerical evidence 
suggests that it is not true in general [n|. Understand- 
ing optimal Hamiltonian simulation and the capacities of 
Hamiltonians in the general case remains an interesting 
open problem. 
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APPENDIX: CATALYTIC SIMULATION OF A 

SUM OF TWO PRODUCT HAMILTONIANS 
WITH THE SAME EXTREMAL EIGENSPACE 

As mentioned in Sec. IIHI Theorem [3] can be extended 
when catalytic entanglement is available. Consider a sum 
of two tensor product Hamiltonians, 

H = J A ® Jb + G a ®G b ■ (89) 

Let J4 denote the restriction of J a to the subspace cor- 
responding to its two extremal eigenvalues. By Rule |3 
J| can be assumed to be traceless. Let G|, Jg, GJ be 
similarly defined. In terms of these Hamiltonians, we 
have 

Corollary 9. Given the resource of catalytic entangle- 
ment, H = J a <8> Jb + Ga <£> Gb is locally equiva- 
lent to [(Aj + Aq)/4] i/ising if the following conditions 
hold: (i) J4 and Gf are supported on the same 2- 
dimensional Hilbert space, and similarly for J B and GJ . 
(it) tr J A e G| =tr J£G§. 

Proof. [(Aj + Ag)/4]iJ Is i ng can simulate H term- wise us- 
ing Lemma |21 and Rule [SI with no need for catalytic en- 
tanglement. 

The following procedure uses H to simulate [(Aj + 

1. Following Rule|Hl Alice and Bob restrict to the ex- 
tremal eigenspace, which is common to both terms 



in H by condition (i). This preserves the extremal 
eigenvalues. The resulting Hamiltonian is essen- 
tially a two-qubit Hamiltonian. 

2. We can assume J4 ® J£ = (Aj/4) ct z a ® o- z b by 
a local change of basis. This can be chosen so 
that G| ® GJ = (A^/4) (cos6>ct z a + sin^cr^) ® 
(cos 9 a z b + sin 0a x s) for some because of condi- 
tion (m). 

3. A further local change of basis takes JJ± ® + 
Gf (g) GJ to its normal form (A 2 / 4) <t Xj 4 ® cr^s + 
(A2/4) cr zA <8> oVb ,6], where A 2 + A 2 = A 2 7 + A 2 G . 

4. Finally, (A|/4) cr^ <8> <J xB + (A 2 /4) o- 2j4 <8> a zB can 
simulate [(A 2 + A^)/4]ct z a ® cr z s using catalytic 
entanglement 10] . 

This completes the proof. □ 

This result allows us to calculate the entanglement ca- 
pacities of the relevant Hamiltonians. If H satisfies con- 
ditions (i) and (ii) of CorollaryEl then 

E H = |(A 2 j + A^). (90) 

There is an input state that achieves Eh without making 
use of ancillas (and in particular, without using catalytic 
entanglement), so Tjj = Eh- As mentioned in Sec. II VI 
this generalizes the case H = \i x o~ x ® a x + [i y o~ y o~ y 
considered in 0. 
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